Numerical study of magnetization processes in rare-earth tetraborides 
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We present a simple model for a description of magnetization processes in rare-earth tetraborides. 
The model is based on the coexistence of two subsystems, and namely, the spin subsystem described 
by the Ising model and the electronic subsystem described by the Falicov-Kimball model on the 
Shastry-Sutherland lattice (SSL). Moreover, both subsystems are coupled by the anisotropic spin- 
dependent interaction of the Ising type. We have found, that the switching on the spin-dependent 
interaction {Jz) between the electron and spin subsystems and taking into account the electron 
hopping on the nearest [t) and next-nearest {t') lattice sites of the SSL leads to a stabilization of new 
magnetization plateaus. In addition, to the Ising magnetization plateau at m,"^ /ml^ — 1/3 we have 
found three new magnetization plateaus located at m,"^ /ml^ = 1/2, 1/5 and 1/7 of the saturated 
spin magnetization m"/ . The ground-states corresponding to magnetization plateaus have the same 
spin structure consisting of parallel antiferromagnetic bands separated by ferromagnetic stripes. 

PACS numbers: 75.10.-b,75.60.Ej,75.40.Mg 
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I. INTRODUCTION 

The Shastry-Sutherland lattice (SSL) was considered 
more than 20 years ago by Shastry and Sutherland [1] 
as an interesting example of a frustrated quantum spin 
system with an exact ground state. It can be described 
as a square lattice with antiferromagnetic couplings J 
between nearest neighbors and additional antiferromag- 
netic couplings J' between next-nearest neighbors in ev- 
ery second square (see Fig. 1). This lattice attracted 
much attention after its experimental realization in the 
SrCu2{BO^)2 compound [1]. The observation of a fasci- 
nating sequence of magnetization (mjma —Xjl, 1/3, 1/4 
and 1/8 of the saturated magnetization m^) in this ma- 
terial stimulated further theoretical and experimental 
studies of the SSL 0,11]. 

As another realization of the SSL the rare-earth tetra- 
borid TmB4 has recently been studied in finite magnetic 
fields [Hi- Since fully polarized state can be reached for 
experimentally accessible magnetic fields, this compound 
allows exploration of its complete magnetization process. 
It was found that the magnetization diagram of TmBi 
consists of magnetization plateaus located at small frac- 
tional values of m/ms =1/7, 1/8, 1/9 ... of the satu- 
rated magnetization, followed by the major magnetiza- 
tion plateau located at m/ms = 1/2. Note that, due to 
large total magnetic moments of the magnetic ions, this 
compound can be considered as a classical system. More- 
over, because of strong crystal field effects, the effective 
spin model for TmB^ has been suggested to be described 
by the spin-1/2 Shastry-Sutherland model under strong 
Ising (or easy-axis) anisotropy Q. From this point of 
view it was natural to begin a description of magnetiza- 
tion process in the TmB^ material from the Ising limit on 
the SSL that can be, in the presence of a finite magnetic 
field h, expressed as follows 



where S'f = ±1/2 denotes the z-component of a spin- 
1/2 degree of freedom on site « of a square lattice and J, 
J' are the antiferromagnetic exchange couplings between 
all nearest neighbor bonds (J) and next-nearest neighbor 
bonds in every second square (J')i as indicated in Fig. 1. 
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FIG. 1: (Color online) The Shastry-Sutherland lattice with 
magnetic couplings J bonds along the edges of the squares 
and J' along the diagonals. 

In spite the relative simplicity of the model Hamilto- 
nian (1), fully different conclusions have been obtained 
for the magnetization curve of this model within vari- 
ous approaches. For example, the authors of Ref. Q 
found, analyzing a finite system consisting of 16 spins, 
a single magnetization plateau at 1/2 of the saturated 
magnetization in accordance with experimental data in 
TmB^. However, numerical simulations obtained within 
the Monte-Carlo and tensor renormalization group meth- 
ods on much larger systems 0, '?] did not confirm this 
conclusion. In contrast to previous results they showed 
that the Ising model on the SSL exhibits in the presence 
of the magnetic field the magnetization plateau only at 
1/3 of the saturated magnetization. Thus the different 
conclusion of Ref. [1] appears to be due to the usage of 
inappropriate finite lattice sizes. 

The existence of the magnetization plateau at only 1/3 
of the saturated magnetization and its absence at 1/2 in- 
dicates that it is necessary to go beyond the classical 
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Ising limit to reach the correct description of the mag- 
netization process in TmB4 and other rare-earth tetra- 
borides. The first such an attempt has been done by 
Meng and Wessel ^ who studied the spin-1/2 easy-axis 
Heisenberg model on the SSL with ferromagnetic trans- 
verse spin exchange using quantum Monte-Carlo and de- 
generate perturbation theory. Besides the magnetization 
plateau at 1/3 of the saturated magnetization they found 
a further plateau at 1/2, which persists only in the quan- 
tum regime. The same results have been obtained by 
Liu and Sachdev analyzing the perturbative effects of the 
transverse fluctuations on the SSL spin multiplets with 
large easy-axis anisotropy Q . 

It should be noted that a similar behavior as for TmBj^ 
has been also observed for other rare-earth tetraborides. 
For example, for ErB4 the rnagnetization plateau has 
been found at m/nis = 1/2 [g, [lO|. for TbB^ at m/nis = 
1/2, 4/9, 1/3, 2/9 and 7/9 [ill and for H0B4 at m/m, = 
1/3,4/9 and 3/5 [To|. 



II. MODEL 

In the current paper we present an alternative model 
of stabilization the magnetization plateaus in the rare- 
earth tetraborides based on the fact that these materi- 
als, in contrast to SrCu2{BO^)2: are metallic. Thus for 
a correct description of ground-state properties of rare- 
earth tetraborides one should take into account both spin 
and electron subsystems as well as the coupling between 
them. Supposing that electron and spin subsystems in- 
teract only via the spin dependent Ising interaction J^, 
the Hamiltonian of the system can be written as 

i 

-hY^{n,^ - n,^) + Hjj, , (2) 



where d^, dicr are the creation and annihilation oper- 
ators of the itinerant electrons in the d-band Wannier 
state at site i and riia = df^di„. The first term of (2) is 
the kinetic energy corresponding to quantum-mechanical 
hopping of the itinerant d electrons between sites i and j. 
These intersite hopping transitions are described by the 
matrix elements tij , which are —t if i and j are the nearest 
neighbors, —t' if i and j are the next-nearest neighbors 
from the SSL and zero otherwise. The second term rep- 
resents the above mentioned anisotropic, spin-dependent 
local interaction of the Ising type between the localized 
spins and itinerant electrons. The third term describes 
an action of the magnetic field on the itinerant electrons. 

To examine the magnetization curve corresponding 
to the model Hamiltonian (2), we have used the well- 
controlled numerical method that we have elaborated re- 
cently to study the ground states of the spinless/spin- 
one-half Falicov-Kimball model [l^ . This method is de- 
scribed in detail in our previous papers [isl . [l^ | and thus 



we summarize here only the main steps of the algorithm: 

(i) Chose a trial spin configuration s = {SI, 5*1, . . . , S'£}. 

(ii) Having s, J^, t and t' fixed, find all eigenvalues 
of h„[s) — tij — aJzSiSij. (in) For a given N — + 
(where N is the total number of electrons) determine 



the ground-state energy E{s) = J^a J2 
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N^) + Hjji of a particular spin configuration s by filling 
in the lowest , one-electron levels . (iv) Gener- 
ate a new configuration s' by fiipping a randomly chosen 
spin, (v) Calculate the ground-state energy E{s'). If 
E{s') < E{s) the new configuration is accepted, other- 
wise s' is rejected. Then the steps (ii)-(v) are repeated 
until the convergence (for given parameters of the model) 
is reached. 



III. RESULTS AND DISCUSSION 

Using the method discussed above we have performed 
exhaustive numerical studies of the model (2) for a wide 
range of model parameters h,Jz,t,t' and J/J' — 1 se- 
lected on the base of experimental measurements [5] . To 
exclude the finite size effects the numerical calculations 
have been done for several different Shastry-Sutherland 
clusters consisting of L = 8 x 8, 10 x 10 and 12 x 12 
sites. The most important result obtained from these 
calculations is that the ground-state spin arrangements 
exhibit the same structure for all examined finite clus- 
ters. In general, this structure is formed by parallel anti- 
ferromagnetic bands separated by ferromagnetic stripes 
and does not depend on the anisotropic spin-dependent 
interaction Jz as well as nearest and next-nearest neigh- 
bor hopping integrals t and t' . The complete list of the 
ground-state spin arrangements (for < m^^ /ml^ < 1) 
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FIG. 2: The complete list of the ground-state spin configura- 
tions (for < m"'' /ml^ < 1) that are stable on finite intervals 
of /i for L = 8 X 8, L = 10 x 10 and L = 12 x 12. The big 
(small) dots correspond to the up (down) spin orientation. 
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that are stable on finite intervals of magnetic field values 
are depicted on Fig. 2. The second important observation 
is that the width w of the antiferromagnetic bands can- 
not be arbitrary, but fulfill severe restrictions. Indeed, we 
have found that with exception the case m'^P/mf — 1/2, 
in all remaining cases the permitted width of the anti- 
ferromagnetic band is only w or w + 2, where w is the 
even number. This fact is very important from the nu- 
merical point of view since it allows us to perform the 
numerical calculations on much larger clusters with the 
extrapolated set of configurations of the above described 
type. The resulting magnetization curves obtained on 
the extrapolated set of ground-state spin configurations 
consisting of parallel antiferromagnetic bands of width 
w {w and w + 2) separated by ferromagnetic stripes are 
shown in Figs. 3-5 for selected values of model parame- 
ters, that represent the typical behavior of the model. 
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FIG. 3: (Color online) Magnetization curves for J'/J — 1, 
Jz = 2, t = 4, = and different values of L. Inset: magneti- 
zation curves of spin and electron subsystems (a = el or sp). 
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FIG. 4: (Color online) Magnetization curves for J'/J — 1, 
Jz = i, t = 4, t' — and different values of L. Inset: magne- 
tization curves of spin and electron subsystems. 

One can see that the switching on the spin-dependent 
interaction between the electron and spin subsystems 
and taking into account the electron hopping on the near- 
est (t) and next-nearest {t') lattice sites of the SSL leads 



to a stabilization of new magnetization plateaus. In ad- 
dition to the Ising magnetization plateau at m^'P/mlP = 
1/3 we have found two new magnetization plateaus lo- 
cated at m^P/mlP = 1/2 and m'P/m'P = 1/5. The 
ground-state spin arrangements corresponding to these 
magnetization plateaus have the same structure consist- 
ing of parallel antiferromagnetic bands of a width w 
(where w = 1 for ra^P /mlP = 1/2, w = 2 for m^P/mf = 
1/3 and w = 4 for m^P /mlP — 1/5) separated by fer- 
romagnetic stripes. Thus, our numerical results show 
that besides the pure spin mechanism (e.g., the easy- 
axis Heisenberg model on the SSL !(^| ) of stabilization the 
magnetization plateaus in rare-earth tetraborides, there 
exists also an alternative mechanism based on the coexis- 
tence of electron and spin subsystems that are present in 
these materials. From this point of view it is interesting 
to compare in detail the ground states obtained within 
these two different approaches. For m^P/mf = 1/3 
our results are identical with ones obtained within the 
Ising d, 0] as well as easy-axis Heisenberg Q model 
on the SSL. The accordance between our and the easy- 
axis Heisenberg solution Q is found surprisingly also for 
/mf = 1/2. In this case both approaches predict the 
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sequence of parallel antiferromagnetic and ferromagnetic 
stripes. For m^P /mf = 1/5 our results postulate a new 
type of spin ordering. 
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FIG. 5: (Color online) Magnetization curves for J'/J = 1, 
Jz = 4, f = 4, f' = QAt and different values of L. Inset: 
magnetization curves of spin and electron subsystems. 

While the magnetization plateaus at m^P/mlP = 1/2 
and 1/3 have been really found in the rare-earth tetra- 
borides d, l9l-[Tl|. the 1/5-magnetization plateau in 
these compounds absent. Instead the 1/5-magnetization 
plateau there have been observed magnetization plateaus 
at smaller values oimfP/mlP, and namely, at m^P/mf = 
1/7, 1/9 and 1/11 [TmBu [5]). Since the sizes of selected 
clusters (60 x 60 and 120 x 120) are not dividable by 
7, 9 and 11 the absence of magnetization plateaus at 
1/7, 1/9 and 1/11 is nothing surprising. To verify the 
possibilities of existence the magnetization plateaus at 
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/mf = 1/7, 1/9 and 1/11 one has to examine much 
larger lattices. Unfortunately, due to the numerical limi- 
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tations we are able to study clusters only slightly higher 
than 120 sites. Such cluster sizes (e.g., L — 140 x 140) are 
sufficient for investigation the stability oim'^P /mf = 1/7 
magnetization plateau, but they are too small for veri- 
fication the magnetization plateaus at ■m''P/mf — 1/9 
and 1/11. In Fig. 6 we present magnetization curves ob- 
tained on clusters consisting of 70 x 70 and 140 x 140 sites 
together with the magnetization curve for L = 120 x 120. 
Comparing these results one can see that a new magneti- 
zation plateau at m'^P/mf = 1/7 is formed and that the 
region of its stability is practically independent of L. 
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FIG. 6: Magnetization curves for J' / J = 1, i = 4, t' = 0, 
Jz = 4 and different values of L. Inset: a detail of magneti- 
zation curves for small h. 

Although we have considered the clusters of different 
classes (the 70 x 70 and 140 x 140 clusters are not di- 
vidable by 3 and the 120 x 120 cluster is not dividable by 
7) the convergence of numerical results to 1/2, 1/3, 1/5 
and 1/7 plateaus is apparent, what indicates that at least 
these plateaus persist in the thermodynamic limit. It is 
not excluded that on much larger clusters also the mag- 
netic plateaus at 1/9 and 1/11 are stable. This conjecture 
supports the fact that the 1/9 and 1/11 phases have the 
same type of the ground-state spin ordering (consisting 
of parallel antiferromagnetic bands (of the even width w) 
separated by ferromagnetic stripes) as the stable phases 
corresponding to the magnetization plateaus at 1/3, 1/5 
and 1/7. However, with respect to our numerical results 
(presented above) showing on the strong suppression of 
the plateau width with decreasing m^^ jml^ we expect 
that the stability region of 1/9 and 1/11 phases will be 
very narrow. 



The magnetization process of the electron subsystem 
is very similar to one described above for the spin subsys- 
tem, but only in the limit m'^P/mf < 0.5 (see Insets in 
Figs. 3-5). Indeed, we have found that for m'^P/mf < 0.5 
the magnetization curves of electron and spin subsystems 
fully coincide for the strong coupling (Jj = 4) between 
electron and spin subsystems, and small deviations are 
observed only for the intermediate coupling (J^ = 2). 
However, a different picture of magnetization processes 
of electron and spin subsystems is observed in the limit 
m^P /m^P > 0.5. In this limit the spin subsystem is al- 
ready fully saturated while the magnetization of the elec- 
tron subsystem changes continuously from m'^'' /mf — 0.5 
to rrf^ /mf = 1. 

In summary, we have presented an alternative model 
of stabilization the magnetization plateaus in rare-earth 
tetraborides based on the coexistence of spin and electron 
subsystems (coupled by the anisotropic spin-dependent 
interaction of the Ising type) in these materials. It was 
shown that the switching on the spin-dependent inter- 
action between the electron and spin subsystems and 
taking into account the electron hopping on the nearest 
and next-nearest lattice sites of the SSL leads to a sta- 
bilization of magnetization plateaus at m^P/m^^ = 1/2, 
1/3, 1/5 and 1/7 of the saturated spin magnetization. 
The ground states corresponding to these magnetization 
plateaus have the same structure consisting of parallel 
antiferromagnetic bands of width ui = 1, 2, 4 and 6 sep- 
arated by ferromagnetic stripes. These results indicate 
that the electron subsystem and its interaction with the 
spin subsystem can play the crucial role in the correct de- 
scription of magnetization processes in rare-earth tetra- 
borides. In our future work we plan to generalize this 
simple model by including the long-range interactions (it 
was shown that such interactions suppress the stability of 
the 1/3 phase [lB|) and considering the Heisenberg spins 
instead of the Ising ones. 
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